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The revival of oscillation and maintaining rhythmicity in a network of coupled oscillators offer an 
open challenge to researchers as the cessation of oscillation often leads to a fatal system degradation 
and an irrecoverable malfunctioning in many physical, biological and physiological systems. Re¬ 
cently a general technique of restoration of rhythmicity in diffusively coupled networks of nonlinear 
oscillators has been proposed in [Zou et ah Nature Commun. 6:7709, 2015], where it is shown 
that a proper feedback parameter that controls the rate of diffusion can effectively revive oscillation 
from an oscillation suppressed state. In this paper we show that the mean-field diffusive coupling, 
which can suppress oscillation even in a network of identical oscillators, can be modified in order to 
revoke the cessation of oscillation induced by it. Using a rigorous bifurcation analysis we show that, 
unlike other diffusive coupling schemes, here one has two eontrol parameters, namely the density of 
the mean-field and the feedbaek parameter that can be controlled to revive oscillation from a death 
state. We demonstrate that an appropriate choice of density of the mean-field is capable of inducing 
rhythmicity even in the presence of complete diffusion, which is an unique feature of this mean-field 
coupling that is not available in other coupling schemes. Finally, we report the first experimental 
observation of revival of oscillation from the mean-field-induced oscillation suppression state that 
supports our theoretical results. 

PACS numbers: 05.45.Xt 


I. INTRODUCTION 

The suppression of oscillation and the revival of oscil¬ 
lation are two opposite but interrelated important emer¬ 
gent phenomena in coupled oscillators. In the suppres¬ 
sion of oscillation, oscillators arrive at a common homo¬ 
geneous steady state (HSS) or different branches of an 
inhomogeneous steady state (IHSS) under some proper 
parametric and coupling conditions. The former is called 
the amplitude death (AD) state [l|, whereas the latter 
is denoted as the oscillation death (OD) state [2[. On 
the other hand, the revival of oscillation is a process by 
which the rhythmic behavior (or rhythmicity) of individ¬ 
ual nodes in a network of coupled oscillators is restored 
from an AD/OD state without changing the intrinsic pa¬ 
rameters associated with the individual nodes Q . 

In the oscillation quenched (suppressed) state the dy¬ 
namic nature of individual coupled oscillators is lost. 
This has many potential applications: For example, the 
AD state is important to suppress unwanted oscillations 
that hinder a certain process, e.g., in laser ^sterns, chat¬ 
tering in mechanical drilling process, etc [l|. Similarly, 
the OD state has a significant role in understanding of 
marw biological processes, e.g., synthetic genetic oscilla¬ 
tor 0 , cardiovascular phenomena [5| , cellular differentia- 
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tion [sl, etc. On the other hand, the research in the topic 
of revival of oscillation is important because many physi¬ 
cal, environmental, biological and social processes require 
stable and robust oscillations for their proper function¬ 
ing: Examples include, El Niho/Southern Oscillation in 
ocean and atmosphere Q, brain waves in neuroscience 
[H, electric power generators [9[, cardiopulmonary sinus 
rhythm of pacemaker cells [T^, etc. In these systems the 
suppression of oscillation may result in a fatal system 
breakdown or an irrecoverable physiological malfunction. 
Thus, it has to be ensured that, if these type of systems 
are trapped in an oscillation cessation state, that state 
has to be revoked in order to establish rhythmicity. 

A recent burst of publications reveal many aspects of 
AD and OD: In particular, identifying coupling schemes 
to induce them [111, [HI , transition from AD to OD 0 , 
their experimental verifications [TsI, etc. But only a few 
techniques have been reported to revoke AD/OD and 
induce rhythmicity in a network of oscillators I, HI- 
In Ref.0 several variants of time delay techniques are 
discussed in order to revoke death states, whereas in 
Ref. 0, network connections are chosen properly in or¬ 
der to revive oscillations. However, most of these tech¬ 
niques lack the generality to revive oscillations from a 
death state. 

Only recently a general technique to revive oscillation 
from the oscillation suppressed state (or death state) has 
been reported by Zou et al. Cl- The authors proposed 
a simple but effective way to revoke the death state and 
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induce rhythmicity by introducing a simple feedback fac¬ 
tor in the diffusive coupling. They showed that this tech¬ 
nique is robust enough to induce rhythmicity in a diffu¬ 
sively coupled network of nonlinear oscillators, such as, 
the Stuart-Landau oscillator, Brusselators model, chaotic 
Lorenz system, cell membrane model, etc. They further 
tested the effectiveness of their proposed technique in 
conjugate M and dynamic coupling [I 3 . 

However, in the absence of parameter mismatch or cou¬ 
pling time-delay, simple diffusive coupling cannot induce 
AD in coupled oscillators. Therefore, for identical os¬ 
cillators under diffusive coupling (without coupling de¬ 
lay) no AD is possible and thus the issue of revoking the 
AD state does not arise. Also, OD in diffusively coupled 
identical oscillators is always accompanied by a limit cy¬ 
cle, thus one needs to choose proper initial conditions 
to revoke that death state. Regarding the conjugate cou¬ 
pling, it is not always a general technique to induce death: 
For example, in a first-order intrinsic time-delay system 
no conjugate coupling is possible. Further, the dynamic 
coupling has its own pitfalls (e.g., its success strongly de¬ 
pends on the intrinsic properties of the oscillators under 
consideration), which has been discussed in detail in [TsI . 

In this context the mean-field diffusive coupling is a 
general way to induce AD/OD even in networks of iden¬ 
tical coupled oscillators: It works in any network of os¬ 
cillators including chaotic first-order intrinsic time-delay 
systems 0- Further, it has been shown in ini that, 
unlike diffusive coupling, the OD state induced by the 
mean-field diffusive coupling is not accompanied by a 
limit cycle. Also, the mean-field diffusive coupling is 
the most natural coupling scheme that occurs in physics 
[m, 0, HqI , biology [3,[il,l2il,[22|, ecology [ 2 ^, etc. Thus, 
for those systems that always need a robust limit cycle 
for their proper functioning, the mean-field diffusive cou¬ 
pling is a much stronger “trap” to induce death in com¬ 
parison with the other coupling schemes. Therefore, it 
is important to study the process of revoking the oscilla¬ 
tion suppression state induced by the mean-field diffusive 
coupling and revive oscillation from the death state. 

Motivated by the above facts, in this paper we intro¬ 
duce a feedback factor in the mean-field diffusive coupling 
and examine its effect in a network of coupled oscilla¬ 
tors. We show that the interplay of the feedback factor 
and the density of mean-field coupling can restore rhyth¬ 
micity from a death state even in a network of identical 
coupled oscillators. Thus, unlike Ref. 0 , here we have 
two control parameters that enable us to revoke the death 
state. Using rigorous eigenvalue and bifurcation analyses 
on coupled van der Pol and Stuart-Landau oscillators, 
separately, we show that the region of the death state 
shrinks substantially in the parameter space depending 
upon those two control parameters. We also extend 
our results to a network consisting of a large number of 
mean-field coupled oscillators and show that the revival 
of rhythmicity works in the spatially extended systems, 
also. Further, for the first time, we report an experimen¬ 
tal observation of the revival of oscillation from death 


states induced by the mean-field coupling that qualita¬ 
tively supports our theoretical results. 


II. REVIVAL OF OSCILLATION WITH 
MODIFIED MEAN-FIELD DIFFUSIVE 

COUPLING: THEORY 

A. van der Pol oscillator 

At first we consider a network of N van der Pol (VdP) 
oscillators interacting through a modified mean-field dif¬ 
fusive coupling; the mathematical model of the coupled 
system is given by 

Xi=yi + e [QX - axi) , (la) 

ili = ai{l - x‘^i)yi - Xi. (lb) 

Here i = 1 • • • and X = ^ Sili the mean-field 
of the coupled system. The individual VdP oscillators 
show a near sinusoidal oscillation for smaller a^, and re¬ 
laxation oscillation for larger Oi . The coupling strength is 
given by e; Q is called the mean-field density mrameter 
that determines the density of the mean-field 0,00 
(0 < Q < 1); it actually provides an additional free pa¬ 
rameter that controls the mean-field dynamics: Q 0 
indicates the self-feedback case, whereas Q ^ 1 repre¬ 
sents the maximum mean-field density. The feedback 
term a controls the rate of diffusion (0 <(a<l): a = 1 
represents the maximum feedback and the original mean- 
field diffusive coupling; a = 0 represents the absence of 
a feedback and thus that of diffusion. Any values in be¬ 
tween this limit can be treated as the intermediate dif¬ 
fusion rate and thus represent a modified mean-field dif¬ 
fusive coupling. The origin of a is well discussed in [HI 
where it is speculated that it may arise in the context of 
cell cycle, neural network and synchronization engineer¬ 
ing. 

As the limiting case we take two identical VdP oscil¬ 
lators: ai ,2 = From Eq.([T]) we can see that there 
are the following fixed points: the origin (0,0, 0,0) and 
two more coupling dependent fixed points: (i) (x*, 

—X*, —y^) where ^* = and y* = ±^Je^a‘^ — 

(ii) {x\ y"', xt, yt) where and = 

±^e2 (a-Q)2 -^(^. 

The eigenvalues of the system at the origin are, 

(a - ea) ± y (a -h ea)^ - 4 

^1,2 — - 2 

_ (g - e(a - Q)) ± y (a -h e{a - Q)Y -1 


From the eigenvalue analysis we derive two pitchfork bi¬ 
furcation (PB) points PBl and PB2, which emerge at the 
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FIG. 1. (Color online) Two {N = 2) mean-field coupled van der Pol oscillators [Eqs. ([T])]: Bifurcation diagram with e for (a) 
a — (b) a — 0.8, (c) a — 0.701 at Q = 0.7. Red (gray) line is for the stable fixed point, black lines are for unstable fixed 

points, green solid circle represents amplitude of stable limit cycle that emerges through Hopf bifurcation and blue open circle 
represents that of an unstable limit cycle, (d) Two parameter bifurcation diagram in the e — Q space for a = 1 (e) Shrinking 
of death region in e — Q space for decreasing a. Area bellow each of the curves for a particular a represents the oscillation 
suppression region and above that shows the oscillating zone, (f) The variation of death region is shown in a — Q space using 
two parameter bifurcation (e = 2). (Other parameter: a = 0.4). 


following coupling strengths: 


^PBl = 


^PB2 = 


a{a - Q) 


( 3 a) 

( 3 b) 


The IHSS, (x*, —X*, —^*), emerges at epBi through 

a symmetry breaking pitchfork bifurcation. The other 
nontrivial fixed point (x’*', , y^) comes into existence 

at epB2, which gives rise to an unique nontrivial HSS. 
Further, equating the real part of Ai^2 and A3^4 to zero, 
we get two Hopf bifurcation points at 

^HBi = —, ( 4 a) 

(a 

( 4 b) 

(a - Q) 

From Eqs. m and o we see that no Hopf bifurcation 
of trivial fixed point occurs for a > 1; in that case, only 
pitchfork bifurcations exist. 

The eigenvalues of the system at the nontrivial fixed 
point {x^, y\ Jx\ Jy^)^ where J = ±1 and j = f or * 
are given by: 




- 4 ciJ 

( 5 a) 

Ai,2 — 

2 


, -62^' ± 1 

A7 

- 4 c 2 J 

5 

( 5 b) 

A3 ,4 = 

2 


where, hi^ = ea — a(l — ci^ = 1 + 2 ax^y^ — aea{l — 

x^‘^), b 2 ^ = e{a — Q) — a(l — C 2 ^ = 1 + 2 ax^y^ — 

ae{a — Q ){1 — x^ ). Now, with increasing Q, eHB2 moves 
towards eppi, and at a critical Q value, say Q*, HB2 
collides with FBI: Q* = a{l — a^). For Q > Q*, the 
IHSS becomes stable at epBS through a subcritical Hopf 
bifurcation, where 


^HBS = y 4 ' ■ (6) 

^a{a-Q) 

This is derived from the eigenvalues of the system at 
the the nontrivial fixed point (x*, —x*, —^*). ^hbs 

actually determines the direct transition from OD to a 
limit cycle, i.e. revival of oscillation. 

The second nontrivial fixed point {x^ ^ y^ ^ x^, y^) that 
was created at epp2 becomes stable through a subcritical 
pitchfork bifurcation at epBS- 


{ 2 a - Q) 

£PBS - • 


( 7 ) 


This is derived from the eigenvalues corresponding to {x ^, 
y^, x^, y^). This nontrivial AD state can also be pushed 
back to a very large value of e by choosing a ^ and 
thus this AD state can also be revoked effectively. 

The above eigenvalue analysis is supported by a nu¬ 
merical bifurcation analysis using XPPAUT [^. Fig¬ 
ure [T](a-c) show the single parameter bifurcation diagram 
depending on e for different a for an exemplary value 
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FIG. 2. (Color online) Network of mean-field coupled van der 
Pol oscillators [Eqs. ([1])] with Q = 0.7, a = 0.4, N = 100: 
Spatiotemporal plots showing (a) and (b) AD and its revival, 
respectively, with decreasing a at e = 2. (c) and d) show the 
OD and its revival, respectively, with a decreasing a at e = 5. 
The upper rows [i.e., (a) and (c)] have a = 1 and the lower 
rows [i.e., (b) and (d)] have a = 0.75. The first t = 5000 time 
is excluded and then x^s for the next t = 100 are plotted. 


valid for a larger network also, we consider the more gen¬ 
eral case oi N = 100 mean-field coupled identical van der 
Pol oscillators of Eqs. m {a = 0.4). Figure. [2J a) shows 
the spatiotemporal plot of the network in the global AD 
regime ata = l, (5 = 0.7 and 6 = 2; here all the nodes 
arrive at the zero fixed point. The global AD state is 
revoked and rhythmicity is restored in the network by 
decreasing the value of a; as shown in Fig. [2][b) for an 
exemplary value a = 0.75. Equivalently, for e = 5 and 
rr = 1 ((5 = 0.7 as before) we get an OD state in the net¬ 
work [Fig. Efc)]. It can be seen that the nodes populate 
the upper and lower branches [shown with yellow (light 
gray) and brown (dark gray) colors, respectively] of OD 
in a random manner and generate a multi-cluster OD 
state. Oscillation in this network is revived by decreas¬ 
ing a; Fig. [2][d) shows this for a = 0.75. Note that the 
values for which AD, OD and oscillations are obtained 
agree with that for the N = 2 case shown in Fig. [TJ 


(5 = 0.7 (throughout the numerical simulation we con¬ 
sider a = 0.4). We observe that the oscillation cessa¬ 
tion state (both AD and OD) moves towards right, i.e., 
a stronger coupling strength e for a decreasing a, and 
for a = 0.701 (^ Q) [Fig. [1] (c)] the death state moves 
much further from PBl in the right direction. We also 
verify that ai a = Q no death state occurs (not shown 
in the figure). We further explore the zone of oscillation 
cessation in a two parameter bifurcation diagram in the 
e — Q space. Figure. [1] (d) shows the bifurcation curves 
for a = 1. The HB2 curve determines the transition from 
oscillation to AD for Q < Q*; beyond this limit the HBS 
curve determines the zone of oscillation and thus that of 
the death region. Figure. [T] (e) shows that the death re¬ 
gion shrinks with deereasing a confirming our theoretical 
analysis. Finally, we plot the phase diagram in Q —a pa¬ 
rameter space at 6 = 2 [Fig. [1] (f)]: we find that a higher 
value oi Q or a lower value of a support oscillations. In¬ 
terestingly, even for a = 1 (i.e., a complete mean-field 
diffusion), one can revive rhythmicity by simply increas¬ 
ing the value of Q; thus, this coupling scheme offers two 
control parameters to revive oscillations. 

Finally, we summarize the observations and discuss the 
following important points: (i) HB2 is the inverse Hopf 
bifurcation point where an AD state is revoked and gives 
rise to a stable limit cycle. This point (or curve in a 
two parameter space) determines the revival of oscillation 
below a critical value (5*, which is determined by a. From 
Eq. ([6|) we see that by choosing a closer to Q {a > Q) 
the death zone shrinks substantially. Thus, to revoke a 
death state one has to choose a ^ Q, and a = Q ensures 
that there will be no death state even in the stronger 
coupling strength (whatever strong it may be), (ii) Even 
if complete diffusion is present, i.e., a = 1, one can still 
achieve the revival of oscillation by choosing Q ^ 1. This 
is an unique feature of the mean-field diffusive coupling, 
which is absent in other coupling schemes. 

To show that our analysis of two coupled oscillators are 


B. Stuart-Landau oscillator 


Next, we consider N Stuart-Landau oscillators inter¬ 
acting through a modified mean-field diffusive coupling 
in their real part; the mathematical model of the coupled 
system is given by 

Zi = (1 + iiOi — \Zi'^)Zi e^QZ — alte(Zi)^ , (8) 

with i = 1 • - N; Z = ^ Sill R^{Zi) is the mean-field 
of the coupled system, Zi = Xi -\- jyi. The individ¬ 
ual Stuart-Landau oscillators are of unit amplitude and 
having eigenfrequency uui. As the limiting case we take 
N = 2, and rewrite Eq. © in the Cartesian coordinates: 


Xi = PiXi — uJiHi + e[QX — axi]^ (9a) 

iji = UiXi + Pitji. (9b) 

Here Pi = 1 - Xi^ - Vi^ {i = 1,2), X = 

We set the oscillators to be identical, i.e., 001^2 = 
From Eq. m it is clear that the system has the fol¬ 
lowing fixed points: the trivial fixed point is the origin 
(0,0,0,0), and additionally two more coupling depen¬ 
dent nontrivial fixed points: (i) (x*, —x*, —y*) where 


* __ ujy 


X = 


uj‘^-\-eay 


and ^* = ± 




2ea 


{x^, y^, x^ ^ y^) where x^ = -- 


. (ii) 
and y^ = 




e{a — Q)y^^-\-uj^ 

e(a — Q) — 2uj‘^ + y^ 

2e(a — Q) 

The four eigenvalues of the system at the trivial fixed 


point (0,0, 0,0) are. 


Ai,2 — 1 — 
A 3 ,4 = 1 — 


e{a — Q) ± 


ea ± y/(eap'^^Tj^ 


(10a) 

(10b) 
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are given by: 


Ai,2 


As,4 


—bi^ ± ybi^ — 4ci-^’ 


62'^ ± Y 62'^ — 4 c 2 -^ 


(13a) 

(13b) 


Where, bi^ = 4(x-^^ + y^‘^) — 2 + e((a — Q), ci-^ = uj‘^ — 
+ (1 - - 3ya^)(l - ea + eQ - Sa;^^ - y^^), 

b 2 ^ = 4(a;a^ + - 2 + ea, C 2 ^' = - 4x^^y^^ + (1 - 

— 3ya^)(l — ea — — y^^). We derive the loci of 

the HBS curve as 


FIG. 3. (Color online) Stuart-Landau oscillators, a; = 2: Bi¬ 
furcation diagram with e (a) a = 1 (Q = 0.5), (b) a = 0.5 
(Q = 0.5). (c) Two parameter bifurcation diagram in e — Q 
space for different a. (d) Spatiotemporal plots of iV = 100 
Stuart-Landau oscillators showing OD at a = 1 (upper panel) 
and the revival of oscillation at a = 0.6 (lower panel): Q = 0.5 
and e = 6. The time scale is same as in Fig. [21 


Through an eigenvalue analysis and also a close inspec¬ 
tion of the nontrivial fixed points reveal that two pitch- 
fork bifurcations (PB) occur at: 


1 +0)2 

^PBl = , 

a 

( 11 a) 

1 + 

epB2 = ^ • 

a - Q 

( 11 b) 


—2((3 -f- (a) + 4>y+ cj^(Qf — Q)[3q!. + Q) 

~ (a-Q)(3a + Q) ’ 

(14) 

Bifurcation diagrams of Fig. Ufa) {a = 1) andjU^b) {a = 
0.51) show that with decreasing a the death region moves 
towards a larger coupling strength e. Figure |3l^c) shows 
this in the e — Q space; one can see that a decreasing a 
shrinks the region of death and thus broadens the area of 
the oscillation state in the parameter space. Finally, we 
consider a network of = 100 identical Stuart-Landau 
oscillators with the same coupling scheme [Eq. ([ 8 ])]. We 
find that the revival of oscillation works here as well. 
Figure |3jd) demonstrates this for an exemplary value of 
e = 6 and Q = 0.5: The upper panel shows the OD 
state in the network for a = 1 , whereas from the lower 
panel we see that rhytmicity is restored in the network 
by reducing the value of a to a = 0 . 6 . 


A symmetry breaking pitchfork bifurcation gives birth 
to the IHSS (x*, ^*, —X*, —y*) at epBi- The second 
nontrivial fixed point {x^^ y\ x’*', y"^) emerges at PB2 the 
stabilization of which leads to a nontrivial AD state that 
coexists with OD. 

Next, we get the Hopf bifurcation point by equating 
the real part of A 3,4 and A 1,2 to zero, 

( 12 a) 
( 12 b) 


= 


a 


^HB2 = 


ol-Q' 


From Eqs. (uni and ([H]) it is clear that for cj < 1 no Hopf 
bifurcations (of trivial fixed point) occur, only pitchfork 
bifurcations govern the dynamics in that case. Erom 
Eq. ([12]) it is noticed that for a fixed a, tuBi is constant, 
but eBB 2 depends only upon Q (but is independent of 
cj, where cj > 1). Now with increasing Q, HB2 moves 
towards PBl, and at a critical Q, say Q*, HB2 collides 
with PBl: Q* = 

The eigenvalues of the system at the nontrivial fixed 
point (xA y^, Jx^ ^ Jy^)i where J = ±1 and ji = f or * 


III. EXPERIMENTAL OBSERVATION OF THE 
REVIVAL OF OSCILLATION 


Next, we implement the coupled system of van der 
Pol oscillators given by Eq. o in an electronic circuit 
(Eigure[l|). We use TL074 (quad JEET) op-amps, and 
AD633 analog multiplier ICs. A ±15 v power supply 
is used; resistors (capacitors) have ±5% (±1%) toler¬ 
ance. The unlabeled resistors have the value R = 10 kfl. 
The op-amp AQ is used to generate the mean-field: 
Vq = which is subtracted by ^ 14 i ,2 

using op-amps denoted by A. One can see that deter¬ 
mines the coupling strength, Rq determines the mean- 
field density and R^ controls the feedback parameter a. 
The voltage equation of the circuit can be written as: 




'Re 


2Rq .a 14,- R, 


R 


E- 

i=i 


_ '±V ■ 

R 


CR 


dVyi 

dt 


R 


VI \ V„i 




(15a) 

(15b) 


Here i = 1, 2. Eqs. m is normalized with respect to 
CR, and thus now becomes equivalent to Eq. ([9|) for the 
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FIG. 4. (Color online) Experimental circuit diagram of the 
modified mean-field coupled VdP oscillators. A, A1-A4, and 
AQ are realized with TL074 op-amps. All the unlabeled resis¬ 
tors have value = 10 kQ. C=10 nF, Ra = 250Q, lA = 0.1 
V. Box denoted by “B” are op-amp based buffers; inverters 
are realized with the unity gain inverting op-amps. G sign 
indicates squarer using AD633. 


following normalized parameters: u 

^ ~ i?, ’ ~ 


i? ’ ^ ^ 100i?a ■ 

Xi = T^, and Hi = Vsat is the saturation voltage 

^sat ^sat 

of the op-amp. In the experiment we take = 0.1 v, 
and C = 10 nF; we choose a = 0.4 by taking Ra = 250 ft 
[using a precision potentiometer (POT)]. 

We experimentally observe the revival of oscillation by 
revoking the oscillation cessation states (AD and OD) 
with varying a (i.e., Ra)- At first we consider the case of 
the AD state: For that we choose e = 2 and Q = 0.7 (by 
setting = 5 kft and Rq =3.5 kD, respectively) and 
decrease a from a = 1 to a lower value. Figure [5][a) shows 
experimental snapshots of the AD state in V^i and Vx 2 at 
Ra = 9.7 kft (i.e., a = 0.97) [using DSO, Agilent make, 
DSO-X 2024A, 200 MHz, 2 Gs/s]; in the same figure we 
show the revival of oscillation from the AD state at an 
exemplary value Ra = 7.07 kfl {a = 0.707). Figure [5][b) 
gives the numerical time series at the corresponding a 
values (Fourth-order Runge-Kutta method with step-size 
of 0.01). The numerical results are in good agreement 
with the experimental observations: Also, the dynamical 
behaviors at these parameter values are in accordance 
with Fig. mi). 

Next, we choose an OD state for e = 5 and Q = 0.7 
(i.e.. Re = 2 kH and Rq = 3.5 kfl, respectively): Ex¬ 
perimental snapshots of the OD state at Ra = 9.42 kQ 
{a = 0.942) and the rhythmicity at Ra = 7.31 kft 
{a = 0.731) are shown in Fig. [5](c). The corresponding 
numerical result is given in Fig. [5][d). We see that the 


Re 


a = 


r = t/RC, 
10^5 = 1 , 
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FIG. 5. (Golor online) (a, c) Experimental real time traces of 
Vxi and Vx 2 along with the (b, d) numerical time series plots 
of xi and X 2 . [(a) and (b)] With = 5 kfl (i.e. e = 2) a 
decreasing Ra (o) restores oscillation (LG) from AD: AD at 
Ra = 9.7 kft (a = 0.97), LG at Ra = 7.07 kft (a = 0.707). 
[(c) and (d)] With Re = 2 kQ (i.e. e = 5) a decreasing Ra 
(a) restores oscillation (LG) from OD: OD at Ra = 9.42 kQ 
{a = 0.942), LG at Ra = 7.31 kft (a = 0.731). Others 
parameters are Rq = 3.5 kfl (Q = 0.7) and Ra = 250 ft 
(a = 0.4). y axis: (a) 200 mv/div (c) 100 mv/div; x axis: 380 
/is/div. 


experimental and numerical results are in good agree¬ 
ment and in accordance with Fig.[T][e). Significantly, de¬ 
spite the presence of inherent parameter fluctuations and 
noise, which is inevitable in an experiment, it is impor¬ 
tant that in both the above cases theory and experiment 
are in good qualitative agreement, which proves the ro¬ 
bustness of the coupling scheme in restoring rhythmicity 
in coupled oscillators. 


IV. CONCLUSION 

We have investigated the effect of a feedback parame¬ 
ter, which controls the diffusion rate, on a network of non¬ 
linear oscillators coupled under mean-field diffusion. We 
have shown that, unlike other coupling schemes, here two 
control parameters exist, namely the density of mean- 
field diffusion and the feedback parameter. The interplay 
of these two parameters can revive rhythmicity from any 
oscillation cessation state: In fact by controlling the feed¬ 
back parameter closer to the density of the mean-field one 
can shrink the region of the oscillation cessation state to 
a very narrow zone in parameter space. More interest¬ 
ingly, even in the presence of complete diffusion (i.e., in 
the absence of feedback parameter), the density of the 
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mean-field alone can induce rhyttimicity from a death 
state. Thus, it offers a very robust oscillation revival 
mechanism. We have extended our study to a network 
consists of large number of nodes and shown that the os¬ 
cillation cessation states can be revoked in that case too. 
Finally, we have supported our theoretical results by an 
experiment with electronic van der Pol oscillators and for 
the first time observed the revival of oscillation from the 
mean-field-diffusion-induced death state. 

Since both the density of the mean-field and the feed¬ 
back parameter have strong connections with many real 
biological networks including cell cycle, neural network, 
synthetic genetic oscillators, etc., we believe that the 
present study will broaden our understanding of those 


systems and subsequently this study will shed light on 
the control of oscillation suppression mechanism in sev¬ 
eral biological and engineering systems. 
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